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Low-temperature transport through molecules with effectively negative charging energy U exhibits a charge- 
Kondo effect. We explore this regime analytically by establishing an exact mapping between the negative-?/ 
and the positive-?/ Anderson models, which is suitable for the description of nonequilibrium transport. We 
employ this mapping to demonstrate the intimate relation between nonequilibrium tranport in the spin-Kondo 
and charge-Kondo regimes, and derive analytical expressions for the nonlinear current-voltage chracteristics 
as well as the shot noise in the latter regime. Applying the mapping in the opposite direction, we elucidate 
the finding of super-Poissonian noise in the positive-?/ Anderson model at high temperatures, by relating the 
correlations between spin flips to pair-tunneling processes at negative U. 
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I. INTRODUCTION 

Recently, several theoretical studies have substantiated that 
transport through single molecules may realize the scenario of 
an Anderson model with negative charging energy U ^2(^1^ 
The negative U originates from the coupling of the electronic 
orbital Ed of the molecule to the molecular vibrations, by a 
mechanism originally introduced by Anderson in the context 
of amorphous semiconductors. 7 Assuming that real excitation 
of the molecular vibrations is energetically not allowed, trans- 
port is described by an Anderson model 

H =sn } j d\d a + Un d ^n dl + ^(e a fe - eV a )c\ k<y c ak a 

a aka 
aka 

with an effectively attractive on-site interaction U < 0. Tun- 
neling of electrons between the molecule and the left (right) 
electrode with amplitude < a =L,fc (t a =R,k) is driven by the ap- 
plied bias V = Vl — Vr. Since only energies in the vicinity of 
the Fermi energy are relevant to transport, we assume the hop- 
ping amplitudes t a independent of k, and a linear dispersion 
relation e ak . 

As a result of the attractive interaction, the system favors 
even electronic occupation numbers and may develop a de- 
generacy between two even-number charge states. At high 
temperatures, finite-bias transport in the vicinity of this de- 
generacy point is accomplished by tunneling of electron pairs. 
In Ref. 6, two of us have demonstrated that this pair-tunneling 
regime can be addressed analytically, and manifests itself in 
distinct features in the linear conductance as well as in the 
nonlinear current-voltage characteristic (IV). In the low- 
temperature regime T < Tk, the degeneracy between the 
two charge states induces electronic correlations and leads to 
a charge-Kondo effect. 8 To date, studies of transport in this 
regime have focused on numerical calculations of the linear 
conductance!^^ 

The central goal of the present paper is to develop an an- 
alytical treatment of linear and nonlinear transport in the 



charge-Kondo regime. We achieve this by establishing a 
one-to-one mapping between the negative-?/ and the conven- 
tional (positive-?/) Anderson model, which we term particle- 
hole/left-right (PHLR) tranformation. The particle-hole trans- 
formation forming the first part of this mapping was intro- 
duced by Iche and Zawadowski, 9 and further elucidated by 
Haldane. 10 This mapping converts between spin and charge 
degrees of freedom of both the localized level and the leads. 
For this reason, this mapping by itself is not well suited for a 
description of nonequilibrium transport. We show that a sub- 
sequent left-right transformation eliminates this problem and 
the resulting combined mapping interchanges charge and spin 
degrees of freedom of the localized level only. We demon- 
strate that this allows for a transparent and analytical inves- 
tigation of transport in the charge-Kondo regime based upon 
well-known results 1 1 for the conventional spin-Kondo effect. 

New physical insight can also be gained by exploiting the 
PHLR transformation in the opposite direction, i.e. translat- 
ing from the negative-?/ to the conventional Anderson model. 
For negative-?/, pair-tunneling processes make it natural to 
expect super-Poissonian current noise. The PHLR mapping 
then implies that Fano factors F > 1 must also occur in the 
positive-?/ Anderson model, where this result is surprising at 
first sight. We show that this noise enhancement can be traced 
back to the influence of a local Zeeman field (which is gener- 
ated by the mapping). This Zeeman field induces correlations 
between consecutive spin-flip processes. 

The presentation of our results is organized as follows. In 
Section||l]we introduce the PHLR mapping, which allows us 
to translate between positive and negative U, and serves as the 
central tool in our subsequent considerations. The conversion 
from negative to positive U is exploited in Section^]] where 
we focus on the occurence of super-Poissonian noise at high 
temperatures outside the Kondo regime. The opposite map- 
ping direction, i.e. from positive to negative U, is employed in 
Section HVl to obtain results for transport in the charge-Kondo 
regime. A summary and discussion of our results are given in 
Section [V] Some additional results on shot noise in the high- 
temperature regime of the negative-?/ model, which are not 
necessary for an understanding of the main text, are relegated 
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to an appendix. 

II. MAPPING BETWEEN THE NEGATIVE- U AND 
POSITIVE-C/ ANDERSON MODELS 

The low-temperature physics of the positive-?/ and the 
negative-[/ Anderson models are intimately related: The for- 
mer gives rise to the conventional spin-Kondo effect^ the 
latter exhibits the charge-Kondo effect. 8 The essential in- 
gredients for the conventional Kondo effect are (i) a spin- 
degenerate localized orbital, and (ii) the SU(2) spin sym- 
metry (including spin-independent hopping amplitudes). By 
close analogy, the central requirements for the charge-Kondo 
effect are (i') degeneracy of two charge states which trans- 
form into each other under particle-hole transformation, and 
(ii') particle-hole symmetry of the entire system including the 
leads. 

This close similarity suggests that a transformation inter- 
changing spin and charge degrees of freedom and thus con- 
verting between SU (2) spin symmetry and particle-hole sym- 
metry may serve as a tool for treating both Kondo effects on 
the same footing. Indeed, as first established by Iche and 
Zawadowski 9 , this is accomplished by a particle-hole transfor- 
mation restricted to one spin direction which maps the charg- 
ing energy U into —U. 

The interchange of spin and charge degrees of freedom has 
several immediate consequences. Tuning the gate voltage in 
the negative-C/ model away from the charge degeneracy point 
causes particle-hole symmetry breaking. This breaking of 
particle-hole symmetry is mapped into a breaking of SU(2) 
spin symmetry in the corresponding positive-?/ model. In- 
deed, the detuning of the gate voltage from the charge de- 
generacy point at negative U maps into a local Zeeman field 
acting on the localized orbital at positive U. Thus, unlike the 
spin-Kondo effect which persists over a wide range of gate 
voltages, the charge Kondo effect fully develops only exactly 
at the charge degeneracy point. Similarly, the spin symmetry 
of the negative-!/ model is mapped into particle-hole symme- 
try of the corresponding positive-C/ model. In the absence 
of a real Zeeman field in the original negative-!/ model, the 
positive-C/ model is thus fixed to the particle-hole symmetric 
point e' d = —U'/2. (Here, we denote quantities belonging to 
the positive-C/ model by primes.) More generally, a real Zee- 
man field only leads to breaking of the particle-hole symmetry 
after the mapping into a positive-C/ model. The charge-Kondo 
effect is thus much less sensitive to (real) Zeeman fields than 
the spin-Kondo effect. 

Accordingly, the central results of the mapping can be sum- 
marized by the following "dictionary", 

e' d = (U + B)/2, U' = -U, B' = 2e d + U, (2) 

where B' denotes the local Zeeman field in the positive-C/ 
model. In extension to Eq. Q, we have also included a pos- 
sible Zeeman field B in the negative-C/ model to emphasize 
the complete formal analogy between the two models. How- 
ever, in the following we will restrict our considerations to 
zero Zeeman field in the original negative-C/ model. We now 



turn to a more detailed derivation of these results, following 
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Particle-hole transformation 

The particle-hole (PH) transformation is carried out for one 
of the two spin directions, and without loss of generality we 
may choose a =\. We define new annihilation and creation 
operators /3, ft by 

[3d] = d|, flak] = Cofcf, (3) 

0di = -d\, aH = c i ( _ fc)r ( 4 ) 

It is straightforward to verify that the (3 operators obey the 
usual anticommutation rules. We define the corresponding 
number operators by n da = ft da (3da and n aka = ft ak<y l3 ak<7 . 
From this, one immediately infers that 

nd] = ndf, n afcT = n afeT , (5) 

ndi = l-n d i, n a ki = 1 - n (_ fc )|, (6) 

and, in particular 

Un d ]n dL = Un d ] - Un d ]n dl . (7) 

The latter equation explains the origin of the sign change of 
the charging energy, and the emergence of a Zeeman field. 
Altogether, the PH transformation maps the original Hamilto- 
nian (up to an irrelevant additive constant) to 

H = y X! " d<T ~ U ^T^rfl + X! ( £afc _ aaeV/2) n aka 

o aha 

+ (t^lkJda + h.c.) + (e d + 17/2) (fidt -n dl ). (8) 

aka 

Here and in the following, we assume symmetric voltage split- 
ting between the leads, i.e. V a = aV/2. [For easier notation, 
we identify a = L,B, with a = +1,-1.] We have also 
chosen to measure k with respect to the Fermi momentum 
kp, so that e a (_A;) = — e ak - Analyzing the resulting expres- 
sion for the Hamiltonian, one notes that the Anderson model 
with negative U has been mapped into a model with positive 
U 1 = —U, which includes an additional Zeeman field B' . 
However, we emphasize that the current through the system, 
originally given by 




now takes the form 




The latter equation clearly demonstrates the crucial problem 
of the PH transformation when considering nonequilibrium 
transport: By the PH mapping, the charge current of the 
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negative-[/ model [Eq. (|9jl] is turned into a spin current in 
the conventional Anderson model [Eq. jlOH. Moreover, fi- 
nite voltages translate into Zeeman field gradients in the trans- 
formed Hamiltonian, see Eq. (|8jl. As a consequence, the PH 
transformation by itself is not a well-suited tool for the analy- 
sis of nonlinear transport. 



Left-right transformation 

A solution is provided by a left-right (LR) transformation, 
performed in addition to the PH transformation. The LR trans- 
formation interchanges the roles of left and right leads for one 
spin direction. In the following, we establish this transfor- 
mation, and show that for symmetric devices (i.e., symmetric 
molecule-lead coupling and symmetric voltage splitting) the 
combined PHLR transformation meets all requirements for a 
one-to-one mapping suitable for nonlinear transport. 

For later convenience, we first carry out the LR transfor- 
mation for general molecule-lead couplings, and only subse- 
quently specialize to symmetric devices. The transformation 
consists of an interchange of left and right labels for the spin 
component previously affected by the particle-hole transfor- 
mation, i.e., 



d' a = fida, 
c Lk.l = 0RH, C R ki 

This transforms the Hamiltonian H into 

U 



@Lkl- 



(ID 

(12) 



aka 



+ 



aka 



( f' r'^ d' 
\ 1 'Qcr aka u a 



h.c 



+ (e d + U/2)(n' d] -n' dl ), 



where 



a] 



t„ 



Rl 



(13) 



(14) 



As noted above, primes denote parameters of the positive-!/ 
model, after the combined PHLR-transformation. The crucial 
point of Eq. dl 31 is that the bias now affects both spin compo- 
nents in the same way, and the current is expressed as 



w 2\dt^ 

\ aka 



(15) 



Through the combined PHLR mapping, the charge current of 
the negative-t/ model [Eq. l|9}] is thus mapped into a charge 
current for the positive-?/ model. 



A. Asymmetric devices in the linear-response regime 

A difficulty arising for asymmetric molecule-lead cou- 
plings consists of the spin dependence of hopping within the 
positive-C/ model, see Eq. (114-1 . In this section, we first show 



that this complication drops out when restricting to the linear- 
response regime. 31 Specifically, we show that the conductance 
of the spin-dependent model H', Eq. dl 31 . is identical to the 
conductance of an ordinary Anderson model 



U 



H " = E n 'd« ~ Un 'd] n 'di + J2( e ak - aeV/2)n' aka 

& aka 

+ E (> c 'ifcX + h - c -) + ( e <* + WKn " n'dO> < 16 > 



aka 



obtained from H' by substituting Eq. J14i with the spin- 
independent prescription t' arJ = t a . 

The calculation of the conductance is simplified by in- 
voking a canonical transformation of the lead operators, 
henceforth referred to as Glazman-Raikh transformation for 
definiteness. 13 Conventionally, i.e. for spin-independent tun- 
neling described by H", this transformation converts the left 
and right lead electrons into two channels ijjika (i = 1,2) ac- 
cording to 





(17) 



where (c, s) 



_ jt R ,t L ) 



-= . This decouples the channel i = 2 

from the localized level, and the tunneling term is transformed 
into 



y^jtipikadq + h.c], 



(18) 



ka 



where t = y/J2a M- As a consequence, the current operator 
can be expressed in terms of the decoupled channel only, 



E | ^tR^lka 1 



h.c. 



(19) 



Turning to the case of spin-dependent hopping described by 
H', we find that the extension of the Glazman-Raikh transfor- 
mation is straightforward, and merely amounts to the intro- 
duction of spin dependences t a — > t aa , t — > t a . The crucial 
point in our proof is that for the spin-dependent tunneling ma- 
trix elements given by Eq. dl4> . the expressions i„ = ^aa 
and tio-t/fcr in fact remain spin independent. As a result, 
the Glazman-Raikh transformation of H' and H" leads to 
the same tunneling terms [Eq. dl SH and current operators 
[Eq. I19H . respectively. Employing the Kubo formalism, 14 
the conductance is expressed in terms of the zero-frequency 
current-current correlation function in equilibrium. Since the 
Glazman-Raikh transformation of H' and H" leads to the 
same expressions for the Hamiltonian and current operators, 
the conductances are identical. 

We emphasize that these arguments do not apply to the 
nonequilibrium case. For finite bias, the different chemical 
potentials of the left and right lead prohibit the channel de- 
coupling via the Glazman-Raikh transformation. 
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TABLE I: Dictionary for the PHLR mapping between the negative-!/ model and the conventional Anderson model with additional Zeeman 
field. Symbols with (without) primes denote quantities in the postive-(7 (negative-?/) model. 




FIG. 1 : Illustration of the mapping between the negative-!/ model (a) 
and the conventional Anderson model (b) with additional local Zee- 
man field B' . All level energies are given as energies per particle. 15 
The horizontal dashed line marks the zero-bias Fermi energy of the 
leads. The one-particle energy Ed is required for adding a single elec- 
tron to the neutral molecule. 




FIG. 2: Fano factor as a function of bias and gate voltage for 
a symmetric junction with Tl = T_r = fcsT, and with U = 
—2 x 10 3 fcsT. Super-Poissonian noise reflects the bunching of elec- 
trons in pair-tunneling processes. 



B. Symmetric devices out of equilibrium 

The spin dependence of the tunneling matrix elements 
Eq. d!4> drops out for symmetric couplings, i.e. when — 
tn we recover an ordinary Anderson model after the PHLR 
transformation. As a result, for symmetric devices we have 
established a one-to-one mapping between transport in the 
negative-?/ model and transport in a conventional Anderson 
model, which remains valid beyond the linear response. The 
conventional Anderson model contains an additional local 
Zeeman field, whose magnitude is determined by the gate- 
voltage detuning from the degeneracy point in the negative-?/ 
model. In addition, spin symmetry of the negative-?/ model 
enforces particle-hole symmetry of the positive-?/ model, 
which is thereby fixed to the symmetric point 2e' d + U' = 0. 
An illustration of the resulting configuration is provided in 
Fig. ^ T ne details of the mapping are summarized by the 
"dictionary" in TableU 



III. FROM NEGATIVE TO POSITIVE U: 
CONSEQUENCES OF A LOCAL ZEEMAN FIELD 

While the principal goal of this paper consists of applying 
the PHLR mapping to study transport in the negative-?/ model 



at low temperatures, it is interesting to note that valuable in- 
sight can also be gained by employing the mapping in the op- 
posite direction. Specifically, in this section we investigate 
the current noise for the positive-?/ model with local Zeeman 
field at high temperatures T Tj^wL 

Within the negative-?/ model, transport in the high- 
temperature regime is accomplished by tunneling of electron 
pairs. 6 Thus, it is natural to expect super-Poissonian zero- 
frequency noise S noise due to electron pairing. A conve- 
nient measure of the effectively transferred charge is given by 
the Fano factor F = 5/2e|(/)|, as demonstrated for Cooper 
pairs in superconductors, 1617 and quasiparticles with frac- 
tional charge in the quantum Hall regimeiAMiSi In Fig. [2] we 
present our numerical results for the negative-?/ model, ob- 
tained via a generalized version of the technique developed 
by Korotkov. 21 (This formalism is dicussed in Ref. |22[ and 
the computation for the pair-tunneling regime is detailed in 
Appendix |X]) 33 The plot depicts the Fano factor of a sym- 
metric device at negative U, as a function of gate and bias 
voltage. As expected, the transfer of electrons in pairs causes 
super-Poissonian noise F > 1 in the regime dominated by pair 
tunneling. The fact that the value F = 2, naively expected 
for electron pairs, is not reached, is explained by the coexis- 
tence various pair-tunneling processes as well as cotunneling 
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FIG. 3: (Color online) Relation between (a) the unidirectional pair- 
tunneling process for negative U and (b) the inelastic spin-flip in the 
conventional Anderson model with local Zeeman field, as obtained 
by the PHLR transformation. The phase space available for the spin- 
up electron in the initial state is marked in red (dark gray). 



of single electrons. The Fano factor for a symmetric device 
typically reaches its maximum in the vicinity of the crossover 
between the cotunneling and pair-tunneling regimes. In this 
region, the phase space for pair tunneling onto the molecule 
(~ eV/2 — Ed — U/2) is large compared to the phase space 
for tunneling off the molecule (~ e d + U/2 + eV/2 ps 0), 
or vice versa. As a result, one junction dominates the pair 
transport with a Fano factor close to 2. A more detailed un- 
derstanding of this behavior may be obtained by an analysis 
of the different processes relevant for the charge transfer. This 
is demonstrated in AppendixlAl 

We now turn to the central point of employing the PHLR 
transformation, and map our negative-?/ results to the conven- 
tional Anderson model with local Zeeman field. It is straight- 
forward to obtain the corresponding expressions for the linear 
conductance and the nonlinear IV as a function of the Zee- 
man field. Remarkably, beyond this the PHLR transformation 
leads to the conclusion that the current noise becomes super- 
Poissonian in the conventional Anderson model. This predic- 
tion comes as a surprise, since charge transfer at positive U is 
essentially comprised of one-particle processes. 

The noise enhancement in the positive-C/ model can be 
traced back to correlations induced by the local Zeeman field. 
Specifically, the origin of the super-Poissonian noise lies in 
the translation of the unidirectional pair tunneling [i.e. paired 
electrons are transferred across the same junction] into the lan- 
guage of the positive- U Anderson model, see Fig. [5] Applying 
the PHLR transformation, one finds that the two consecutive 
pair processes correspond to two inelastic cotunneling pro- 
cesses: The first process flips the spin from down to up, the 
second process flips it back to a =J.. In both cases, one elec- 
tron is transferred from the left to the right lead. The crucial 
point is that these two processes are correlated whenever the 
Zeeman field is nonzero. For B' > 0, a spin flip from a ={ to 
a =t requires (twice) the Zeeman energy. The same energy is 
released for the opposite spin flip. As a result, the phase-space 
involved in these inelastic processes causes one spin flip to be 
faster than the other one. The difference of the rates becomes 



maximal when the Zeeman field becomes as large as the bias 
voltage, B' « eV. These correlations cause the enhancement 
of noise beyond the Poissonian limit. 



IV. FROM POSITIVE TO NEGATIVE U: ANALYTICAL 
STUDY OF THE CHARGE-KONDO EFFECT 

In this section, we employ the PHLR mapping to study 
transport in the negative- U model, both in the poor-man- 
scaling regime and the fully developed Kondo regime. 

A. The onset of the charge-Kondo effect: Logarithmic 
corrections 

The degeneracy between the two charge states n — and 
n = 2 in the negative-?/ Anderson model results in Kondo 
correlations at low temperatures. 8 For temperatures T ^> Tk, 
the development of these correlations is expected to be sig- 
nalled by logarithmic corrections to the leading-order pertur- 
bative results. Here, we exploit the PHLR transformation 
to extract the logarithms for the negative- U case from well- 
known results in the context of the conventional Anderson 
model, see e.g. Ref.lTll 

Specifically, we employ the following chain of transforma- 
tions. Our starting point is the negative-?/ Anderson model. 
Applying the PHLR transformation, we obtain a conventional 
Anderson model with local Zeeman field, Eq. Jl 3i . Then, a 
Schrieffer- Wolff transformation is performed. 23 Up to an ir- 
relevant renormalization of the level energy Ed, the mapping 
results in H' = H + Hj + Hz with 

H a = ^(e afc - aeV/2)cl ka c aka , (20) 

aka 

Hj = ^ Jaa'kk'&aa'kk' ' S, (21) 
aa' kk' 

H z =B'S Z . (22) 



Here, the exchange constant J and the lead electron spin op- 
erator s are given by 



T X< J.' 

J aa kk' — "a a 



1 



1 



e a fe 



Ca'fc' — £ r . 



w 



V 



£ d ~ e a >k> 



and 



aa' kk 



( c LcT Ca ' fc 'T _ c L| c «'fc'.l)A 



kk' ~ C ak\ Ca ' k 'ii 



'kk' 



C akl Ca ' k 'T- 



(23) 

(24) 
(25) 



Note that the usual potential scattering term in H' vanishes 
due to the particle-hole symmetry of the positive-?/ model. 
As usual, the exchange is approximated by its Fermi-energy 
value, Jaa'kk' J aa' = — &t a t a > /U within a band of width 
D and zero outside. In the linear-response regime, the Hamil- 
tonian may be further simplified by diagonalizing J = (J aa ') 



6 



with eigenvalues and J = — (8 /U) J2 a \t a \ -In terms of the 
eigenvectors of J, this results in the one-channel problem, 11 



(26) 



Alternatively, this can be derived by applying the Glazman- 
Raikh transformation 13 to the Anderson Hamiltonian before 
carrying out the Schrieffer- Wolff transformation. 

After the mapping of the negative-C/ problem to the Kondo 
Hamiltonian, we are now ready to translate results from 
the conventional spin-Kondo effect to the charge-Kondo sce- 
nario of negative U. In particular, a perturbative treatment 
of the exchange coupling leads to the well-known logarith- 
mic terms ~ J 2 ln(D/ui) in the amplitude for a transition 



ki 



11.12 



Summing up the leading loga- 



rithmic contributions to all orders, or using poor-man's scal- 
ing alternatively, 24 the exchange coupling is renormalized as 
v.] -> l/\n{T/T K ), where T K = Dexp[-l/uJ] denotes 
the Kondo temperature, and v the density of states of the leads. 
Employing the Kubo formula, the conductance for vanishing 
Zeeman field B' — 2ed + U = is found to be 1 1 



G = 



2e 2 4T L r R 3tt 2 



1 



h (Tl+Tr) 2 16 \n 2 (T/T K ) 



(27) 



tion of the Kondo fixed-point Hamiltoniani 2 ^ 2 ^ given by 

ka 

k,k' a 

^T^T^l^i: 



+ 



~kv 2 Tk 



(30) 



where ipka annihilates a quasiparticle with energy £/. = vpk 
and spin a. The a-term describes the quasiparticle scattering 
at the impurity site, the /3-term accounts for the induced inter- 
action between quasiparticles. As shown by Nozieres, the fact 
that the Kondo resonance is floating on the Fermi sea, fixes the 
parameter ratio a/ (3 to unityi2^2& The additional 7-term de- 
scribes the effect of a local Zeeman field B' acting on the im- 
purity (as generated by the PHLR mapping). We note that the 
coefficients a and 7 take on the values a = 7 = 1, when we 
interpret Tk in Eq. | |30I > to be the exact Kondo temperature, as 
determined e.g., from the exact Bethe-ansatz solution 
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i.e., 



T, 



K 



rpba 
1 K 



(2UT) 



1. Linear conductance 



Note that the leading order correctly reproduces the rate- 
equations result for the pair-tunneling peak: For T ^> Tk 
we have \n(T/T K ) ~ 1/vJ, so that 



G w 2Ae 2 T L T R /U 2 h, 



(28) 



in agreement with results in Ref. For decreasing tempera- 
ture, the developing Kondo correlations cause a slow logarith- 
mic increase of the peak height value. 

Additional corrections affect the tails of the conductance 
peak. In the language of the positive-C/ Anderson model, this 
implies considering the situation of a large Zeeman field B' = 
2sd + U ^ ksT . For this limit, a relation similar to Eq. J27I 
can be derived, 1 1 namely 



G = 



2e 2 ATrF 



Li- R 



1 



h {T L + T R ) 2 16 \n 2 [(2e d + U)/T K }' 



(29) 



This captures the logarithmic corrections to the cotunneling 
tails of the conductance peak in the negative-C/ model. 



B. The fully developed charge-Kondo regime 

We now turn to the investigation of the fully developed 
Kondo regime in the low-temperature limit, T -C Tk, and em- 
ploy the PHLR transformation to translate between the con- 
ventional spin-Kondo effect and the charge-Kondo effect. Our 
considerations are based on Nozieres' Fermi-liquid descrip- 



Within the Fermi-liquid description, the linear conductance 
may be expressed as a function of the scattering phase shift. 
As shown by Nozieres in Ref. 25, this phase shift depends 
both on the energy of the incoming particle and on the quasi- 
particle distribution, and can be expanded as 



7T 

2 



ae 
Tk 



(3n s jaB' 
vT K T K 



(31) 



where we denote f =j, J =|. The leading-order phase 
shift 7r/2 makes the molecule a perfectly open channel in the 
Kondo regime. From the phase shift we can determine the 
transmission and reflection amplitudes from left to right (from 
right to left), t„ and r a (t' a and r' a ), respectively. Indeed, us- 
ing the Glazman-Raikh transformation, 13 one obtains the cor- 
responding scattering matrix 



S«r(e) 




This results in a transmission coefficient 




(32) 



4|*l| 2 M 2 
(M 2 + M 2 ) : 



■sin 2 (<5 CT ). 



(33) 



At T = 0, the linear conductance is not expected to depend 
on a or (3, since the quasi-particles have energy e = in lin- 
ear response, and n„ = in the absence of a global magnetic 
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field. As a result, the Kondo-regime conductance of the An- 
derson model with local Zeeman field B' <C Tjc is given by 



2e 2 4r f ,r 



ft (rL + r fl ) : 



[I-^B'/TkY] 



(34) 

Exploiting the identity of linear conductances for the positive 
and negative-?/ models, see Section Hi Al we conclude that 
the charge-Kondo effect in the negative-?/ model leads to the 
conductance 



G = 



2e 2 4T L T! 
h {Tl + T, 



(35) 



Here, the correspondence between Zeeman field (at positive 
U) and gate voltage (at negative U) is directly reflected in 
a departure from the unitary limit as soon as the negative-?/ 
system is tuned away from the charge-degeneracy point. 



2. Nonlinear current-voltage characteristics 

For the special case of symmetric junctions, the PHLR 
transformation allows us to go beyond linear response. For 
nonzero bias, quasiparticles with finite energy are involved, 
and hence the a- and /3-terms in the effective Hamiltonian, 
Eq. (I30i . become relevant. These terms describe the weak 
scattering off the spin-singlet state as well as the induced 
interaction between quasiparticles, and cause a reduction of 
the unitary current I u = 2e 2 /hV due to backscattering 
events. The total current may thus be written in the form 
I = I u — lb- In order to extract the backscattering contri- 
butions from the effective Hamiltonian, we relate the quasi- 
particle states ipka to left-movers and right-movers, L^ a and 
Rka, respectively. For a symmetric device, this relation is 
given by tpka = ^i^ka + Rka)- M Importantly, left-movers 
originate from the right lead and right-movers from the left 
lead. This fact allows one to account for the nonequilibrium 
situation at finite bias by identifying the distribution of right- 
movers (left-movers) with the Fermi distribution of the left 
(right) lead, i.e. /l(_r)(c)- Substituting the LR decomposition 
into the effective Hamiltonian, one finds that the elastic term 
(~ a) and the inelastic term (~ (3) generate backscattering 
transitions which turn left-movers into right-movers, and vice 
versa. Since these contributions act as weak perturbations, one 
may evaluate the backscattering current by summing the cor- 
responding backscattering rates obtained via Fermi's golden 
rule. 

Following Ref. 28 and including the additional local Zee- 
man term results in the total current 



2e 5 



V 



1 - ^B'/Tk) 



5/? 2 



12 



V 



(36) 



Due to the onset of backscattering at finite energies, the cur- 
rent is reduced with increasing bias voltage. The breaking 
of spin-symmetry by the local Zeeman field leads to an ad- 
ditional reduction ~ (B') 2 , as required by symmetry. Now, 



we return to the charge-Kondo effect in the negative-?/ model 
by applying the PHLR transformation. The resulting current 
close to the unitary limit is given by 



2e 2 



V 



1 - ( 7 [2e d + U]/T K f 



12 



5/3 2 ( V_ 



(37) 

revealing the current reduction due to a gate detuning from the 
charge degeneracy point. 



3. Shot noise 

We finally turn to a discussion of shot noise in the Kondo 
regime, which we have recently investigated for the conven- 
tional Anderson model in Ref. 28. The PHLR transformation 
allows us to transfer our central results for the conventional 
Anderson model to the negative-?/ model. 

As demonstrated in the previous subsection, close to T = 
and for small bias voltages, the current is nearly unitary, and 
backscattering events are rare. In this scenario, a sensible def- 
inition of the Fano factor does not involve the transmitted cur- 
rent (which would yield F ~ 0), but rather the backscattering 
current 



F = S/2e\I b 



(38) 



The backscattering current consists of a competition be- 
tween single-particle and pair backscattering processes. Re- 
markably, for vanishing local Zeeman field and the conven- 
tional positive-?/ Anderson model, we have shown that the 
Fano factor is super-Poissonian, F = 5/3, due to the pair 
backscattering. 28 It is the enhanced phase space for pair scat- 
tering events which renders their contribution dominant, de- 
spite the fact that they constitute only one of the seven rele- 
vant processes. 28 The inclusion of an additional local Zeeman 
field is straightforward, and we obtain 



F = 



2{B'/T K f 



W/Tk? 



2{B'/T K y + {V/T K f 



(39) 



(Here, we have exploited the identity of the coefficients a and 
7 following from the Bethe ansatz. 27 ) Hence, for B' = we 
reproduce the result F — 5/3, while in the opposite limit of 
large Zeeman fields B' 3> V, single-particle backscattering at 
the singlet dominates and the Fano factor is Poissonian, F = 
1. 

Applying the PHLR mapping now yields the corresponding 
Fano factor for the negative-?/ model, 



F = 



2([2e d + U]/T K f 



W/Tk? 



2{[2e d + Uj/Z 



K 



K 



(40) 



The Fano-factor enhancement to 5/3 only dominates the im- 
mediate vicinity of the charge-degeneracy point. A detuning 
from this point results in a suppression of the Fano factor, and 
Poissonian noise is recovered as soon as the detuning is large 
compared to the bias voltage. 
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V. SUMMARY 

By an extension of the particle-hole transformation due 
to Iche and Zawadowski, 9 we have established a one-to-one 
mapping between the negative-?/ and the positive-?/ Ander- 
son models appropriate for the situation of nonequilibrium 
transport. This mapping transforms spin degrees to charge de- 
grees of freedom, converts between SU(2) spin symmetry and 
particle-hole symmetry, and is thus ideally suited to the inves- 
tigation of the spin- and charge-Kondo effects. In the case of 
symmetric devices, the mapping leads to a conventional An- 
derson model with spin-independent hopping terms. For this 
scenario, we have demonstrated the usefulness of the mapping 
by deriving analytical expressions for nonlinear transport in 
the negative-?/ model, both in the poor-man-scaling and the 
fully developed Kondo regimes. 

Applying the mapping in the opposite direction, we have 
shown that an additional Zeeman field causes correlations in 
the high-temperature transport through a positive-?/ device. 
These correlations are reflected in super-Poissonian noise, and 
can be directly related to pair-tunneling processes at negative 
U for which enhanced Fano factors occur naturally. 

In closing, we emphasize that the nonequilibrium physics 
in the charge-Kondo regime remains a partially open ques- 
tion for asymmetric molecule-lead couplings. Specifically, 
for asymmetric couplings between the molecule and the leads, 
the hopping terms in the corresponding positive-?/ Anderson 
model become spin dependent, a situation not commonly ad- 
dressed in the Kondo literature to date. 
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APPENDIX A: HIGH-TEMPERATURE NOISE 
CALCULATIONS 

Our calculations of the current shot noise for pair tun- 
neling are carried out by means of a generalized version of 
Korotkov's formalism, 21 which is detailed in Ref. |22j This 
method is based on a decomposition of the time-dependent 
current into discrete contributions from individual tunneling 
processes. The relevant processes labelled by v = 1, . . . , 5 
are identified in Fig.|4] In addition to the transition rates, the 
formalism requires as input the charge sfj v (in units of e), 
transferred across junction a for a process of type v with ini- 
tial and final state % and /, respectively. For completeness, 
these charge-transfer coefficients are summarized in Table HP 



s 02,i/ 


a — L 


a = R 


s 20,i/ 


a = L 


a = R 


V = 1 


2 





V = 1 


-2 





2 





-2 


2 





2 


3 


1 


-1 


3 


-1 


1 


a 

S W.L> 


a — L 


a = R 


~ a 

S22,v 


a = L 


a = R 


4 


1 


1 


4 


1 


1 


5 


-1 


-1 


5 


-1 


-1 



TABLE II: All nonvanishing charge-transfer coefficients sft v char- 
acterizing the current contributions of individual processes. Sketches 
of the various processes u = 1, . . . , 5 are depicted in Fig.0] 



Note that in distinction to Ref. 22, the absolute coefficients 
obtain values larger than unity due to pair tunneling. 




V=4 x2 



v=5 



FIG. 4: Elementary processes for the transitions | ) — > | 2 ) and 
| ) — * | ). The label "x2" signals an additional factor of 2 due 
to two incoherent spin contributions. The diagrams for the processes 
corresponding to the transitions | 2 } — > j ) and | 2 ) — > | 2 ) are 
obtained by reversing the direction of all arrows. 

It is also possible to obtain a comprehensive analytical un- 
derstanding of the numerical results based on the Fano factor 
formula derived in Ref.|29j 



F = (Ni) 



(Al) 



In other words, the Fano factor acquires two distinct contribu- 
tions: The first term in Eq. JAU arises from the fluctuations 
of waiting times t.- L , the second term reflects the fluctuations of 
the number of transferred electrons JV,-. Equation JA1I is valid 
whenever the quantities Ni and t{ are uncorrected, and have 
fixed probability distributions independent of i. Employing 
this formula, we analytically calculate the separate Fano fac- 
tors for cotunneling, unidirectional pair tunneling, and mixed 
pair tunneling in the following. 



1. Cotunneling 

For temperatures small compared to the applied bias, 
ksT <C eV, cotunneling transport is essentially unidirec- 
tional: transitions with electron transfer in the direction op- 
posed to the applied bias can be neglected. The current is then 
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determined by a single parameter, the cotunneling rate W. 
The evaluation of Eq. flAll > is straightforward. The number 
of transferred electrons per event is exactly one, and the fluc- 
tuations of the transmitted charge vanish. The waiting times 
ti follow an exponential distribution p(t) = We~ wt , so that 
the first and second moments are given by 



<*?> = 



w 2 



As a result, cotunneling leads to a Fano factor of unity, 

(g) - (U) 2 



(A2) 



(A3) 



This explains the Poissonian noise in the regions dominated 
by cotunneling in Fig.|2| 



2. Unidirectional pair tunneling 

We now consider the case of unidirectional pair tunneling, 
i.e. transport is assumed to be caused by the following se- 
quence of events: (i) An electron pair tunnels from the left 
lead onto the molecule, (ii) The electron pair tunnels off of 
the molecule into the right electrode. When evaluating the 
Fano factor via Eq. ( I Alt , we consider events in one of the 
two junctions. Without loss of generality, we choose the left 
junction. Then, each event (i) transfers N = 2 electrons, and 
there are no fluctuations of this number. However, it is crucial 
to note that the waiting time ti now consists of a sum of two 
times, the waiting time for the event (i) in junction L and the 
waiting time for event (ii) injunction R, ti = tf + tf. Each 
waiting time tf' R is exponentially distributed according to 



Due to their phase-space behavior, the pair-tunneling rates 
vary with gate voltage according to 

W$L 2 ~ X (et//2 ~ £d ~ U / 2 ) e ( eV / 2 - £ d- V l% 

(A10) 

~^f( £ d + u / 2 + eV/2)6(e d + U/2 + eV/2), 

(All) 

valid in the limit of low temperatures. Consequently, for a 
symmetric device the rates are identical at the degeneracy 
point 2e d +U = 0, and the Fano factor resulting from Eq. dA9l > 
is F = 1. Away from the degeneracy point, one rate loses 
phase space while the other gains phase space. As a result, 
for approximate alignment of the two-particle level with one 
of the Fermi energies, Ed + U/2 w ±eV/2, one rate nearly 
vanishes and the other one remains finite. The corresponding 
Fano factor is F = 2. This dependence of the Fano factor 
on gate voltage is clearly reflected in the numerical results 
depicted in Fig. [2] Interestingly, for asymmetric devices uni- 
directional pair-tunneling leads to Fano factors F > 1 even at 
the degeneracy point. In the limit r a /r a / — * oo, one obtains 
a Fano factor of F — 2. 

It is worth noting that Eq. < IA9I > also yields the Fano factor 
for sequential tunneling when replacing (N) — 2 by (Ni) = 
1, as well as W$U 2 -> W$L>i and W 2 R _^ -> Wft> . The 
crucial difference between sequential and pair tunneling is that 
sequential rates are independent of gate voltage as long as the 
level position remains in the bias window. Accordingly, the 
Fano factor for sequential tunneling in a symmetric device is 
equal to 1/2. 



PL(t)=W t 2 exp[-W^ 2 t] 
PR{t)=W^ exp[-Wi^t] 



(A4) 
(A5) 



where Wq_^ 2 (WzUo) is the rate for the pair-tunneling tran- 
sition in the left (right) junction. The probability distribution 
P(t) for the sum of waiting times in the left and right junc- 
tions is simply given by the convolution 



P(t) 



dt L / dt R p L (t L ) PR (t R )5(t -t L - t R ) 



-Wf 



(A6) 



The resulting first and second moments of the total waiting 
time are 



(t 



1 



w L _ 2 + wH 



wH 



W L W R 



W^ 2 ?[w R ^Y 

Substituting into Eq. (IA1> . we obtain for the Fano factor 



(A7) 
(A8) 



(ti 



(w R ^ + w^ 2 y 



3. Mixed pair tunneling 

Finally, we consider mixed pair tunneling, i.e. the trans- 
port mode typical for pair tunneling in asymmetric junctions. 
Two electrons enter the molecule from the left lead (rate Wi), 
and subsequently split up to leave the molecule via the left 
and right junction, respectively (rate W 2 ). Let us consider the 
right junction. Here, exactly one electron is transferred for 
each combined process of pair tunneling onto and off of the 
molecule. The waiting time again consists of a sum of two 
exponentially distributed random times, exactly as in the uni- 
directional pair-tunneling case. All arguments given there can 
be reapplied for the case of the mixed pair-tunneling case, re- 
sulting in 



F 



(W 1 + W 2 f 



(A12) 



Formally, this is identical to the result for sequential tunnel- 
ing. Once again, the important difference is the gate-voltage 
dependence of the pair-tunneling rates. While the phase-space 
behavior of the rate W\ is given by Eq. JA10> . the splitting 
pair rate scales like 



II, 



r r ,r 



Li R 



(2e d + U)9(2e d + U). 



(A13) 
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processes. Outside the pair-tunneling regime, cotunneling is 
the only relevant process and leads to a Fano factor equal to 
1. Inside the pair-tunneling regime, cotunneling, as well as 
unidirectional and mixed pair tunneling coexist. The resulting 
Fano factor is given by a weighted average of the individual 
Fano factors. This comprehensively explains the qualitative 
features of the numerical results in Fig. [2] and the fact that 
the Fano factor does not fully reach the value of 2 naively 
expected for pair tunneling. 



e d + \U\ /2 (1(T 3 \U\) e d + \U\ /2 (1CT 3 \U\) 



FIG. 5: Fano factor as a function of gate voltage at finite bias 
ieV — 50 x 10 -3 \U\) for different magnitudes of junction asymme- 
try r_L/r_R. The dominance of mixed pair-tunneling processes for 
asymmetric junctions reduces the super-Poissonian noise and leads 
to Fano factors below 1 . 

Consequently, the Fano factor for the mixed pair-tunneling 
process varies between F — 1 and F = 1/2, depending on 
the coupling ratio and gate voltage. 

4. Interpretation of the full Fano factor 

The full Fano factor, Fig. [2] can now easily be interpreted 
in terms of the contributions from the previously discussed 



For completeness, we briefly comment on the consequences 
of asymmetric coupling Tl ^ T/j on the shot noise in the 
negative-C/ model. As discussed for current-voltage charac- 
teristics in Ref. 6, the degree of coupling asymmetry plays 
an important role in determining the relevant transport pro- 
cesses. As a result, the noise characteristics crucially depend 
on the ratio Tr/Tl, see Fig. [5] These plots show results for 
the finite-bias Fano factor as a function of gate voltage, and 
for different coupling ratios. For devices with large asymme- 
try, r a /r a ' ^> 1, unidirectional pair-tunneling is suppressed 
and mixed pair processes take over. 6 As a result, the Fano fac- 
tor is reduced, and sub-Poissonian noise dominates the trans- 
port, see Fig. [5] The influence of unidirectional pair-tunneling 
close to the degeneracy point is reflected in remaining super- 
Poissonian traces for moderate asymmetries. 
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